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ABSTRACT 

The spatial distribution of the Galactic satellite system plays an important role in 
Galactic dynamics and cosmology, where its successful reproduction is a key test of 
simulations of galaxy halo formation. Here, we examine its representative nature by 
conducting an analysis of the 3-dimensional spatial distribution of the M31 subgroup 
of galaxies, the next closest system to our own. We begin by a discussion of distance 
estimates and incompleteness concerns, before revisiting the question of membership 
of the M31 subgroup. We constrain this by consideration of the spatial and kinematic 
properties of the putative satellites. Comparison of the distribution of M31 and Galac- 
tic satellites relative to the galactic disks suggests that the Galactic system is probably 
modestly incomplete at low latitudes by ~ 20 %. We find that the radial distribution 
of satellites around M31 is more extended than the Galactic subgroup; 50 % of the 
Galactic satellites are found within ~ 100 kpc of the Galaxy, compared to ~ 200 kpc 
for M31. We search for "ghostly streams" of satellites around M31, in the same way 
others have done for the Galaxy, and find several, including some which contain many 
of the dwarf spheroidal satellites. The lack of M31-centric kinematic data, however, 
means we are unable to probe whether these streams represent real physical associ- 
ations. Finally, we find that the M31 satellites are asymmetrically distributed with 
respect to our line-of-sight to this object, so that the majority of its satellites are on 
its near side with respect to our line-of-sight. We quantify this result in terms of the 
offset between M31 and the centre of its satellite distribution, and find it to be signif- 
icant at the ~ 3 a level. We discuss possible explanations for this finding, and suggest 
that many of the M31 satellites may have been accreted only relatively recently. Al- 
ternatively, this anisotropy may be related to a similar result recently reported for 
the 2dFGRS, which would imply that the halo of M31 is not yet virialised. Until such 
time as a satisfactory explanation for this finding is presented, however, our results 
warn against treating the M31 subgroup as complete, unbiased and relaxed. 
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1 INTRODUCTION 



The Galactic satellite system consists of approximately 
a dozen galaxies. This is an order of magnitude less 
than the number of haloes pr edicted to exist as sate llites 
from CDM simulations (e.g. iKauffmann et aD Il993l) . al- 
though this conflict is potentially resolved by suppressing 
star formation in s ome subset of dark matter haloes (e.g. 
iBullock et afl feuOOb Successful reproduction of the Galac- 
tic satellite system, in terms of internal properties, fre- 
quency, and spatial properties, is considered an important 
test of cosmological galaxy formation models. In galactic 
dynamics, to o, satellite populations pla y an important role . 
For example. ILitt le fc Tr emainei Jl 987lb iKo chanekl Jl996f) . 
IWilkinson fc Evans! Jl999h . lEvans fc WilkinsonI ll2000h and 



lEvans et alJ <l200 Ql have used the Galactic and M31 satellite 
systems to calculate the dynamical mass of the host galaxies. 

IWillman et al.l (|2004h have questioned whether the cur- 
rent census of Galactic satellites is complete, and suitable for 
use in the above ways. They conclude that several satellites 
are probably missing at low Galactic latitudes, and suggest 
that searches for Galactic satellites could be suffering from 
incompleteness issues at large Galactocentric distances. Ac- 
cording to their study, these effects could have resulted in 
a factor of up to 3 discrepancy between the observed and 
actual number of Galactic satellites. 

Around distant galaxies, a preference for satellites to lie 
along the minor axis of the host galaxy has been reported, 
implying a preference for polar orbits. Th is is th e so-called 
Holmberg ef fect, first noted by jHolmberd <|l969D and later 
revisited bv IZaritskv et all dl997l) . However, recent results 



2 McConnachie & Irwin 



from SDSS question these findings, and i ndicate a prefer - 
ence for orb i ts in the plane of the disk jBrainerdlEool . 
iKnebe et alJ l)2004t) suggest these anisotropics are the result 
of the accretion of satellites along filaments, leading to a 
preference for satellites to be observed along the major axis 
of the host halo and reflecting their initial infall onto the 
host galaxy. However, given the ambiguous and contradic- 
tory nature of the observational studies, the significance of 
this is not clear. 

The Holmberg effect is not observed around nearby sys- 
tems, but the Galactic satellite distribution is anisotropic. 
The LMC and SMC are on the opposite side of the sky to 
the Ursa Minor and Draco dwarf spheroidals (dSphs) , on an 
almost polar great circle on the sky. On the same gre at circle 
lies the Magellanic Stream iMathewson et al"lll974l) . a large 
trail of HI gas which has been stripped from the M agellanic 
Clouds. This correlation in position led iLvnden-B cll (1976) 
to hypothesise that the LMC, SMC, Ursa Minor and Draco 
were once part of a "Greater Magellanic Cloud" , which was 
torn apart in the gravitational potential of the Galaxy and 
condensed to form the latter two galaxies. The LMC and 
SMC are proposed to be the sur viving remnant of t his body. 

At around the same time, iLynden-Bell Jl982l) noticed 
that the Fornax, Leo I, II and Sculptor dwarf galaxies 
all appeared to align in another polar great circle (the 
"FLS" stream), in a similar fashion t o the Magellanic 
Stream. Later work by Maicwski il994T) showed that the 
then newly discovered Sextans and Phoenix dwarf galaxies 
both aligned with this plane, along with most of the Galac- 
tic globular clusters with the red dest horizontal branches. 
iLvnden-Bell fc Lvnden-Belll il995T> searched for all possible 
streams of satellites in the halo of the Galaxy, including some 
of the distant globular clusters, and found several candidates 
in addition to the Magellanic and FLS streams. 

Kinematic evidence for these streams is still inconclu- 
sive. Only radial velocities can be measured readily and 
although some proper motions for these satellites exist, a 
common cons ensus has yet to em erge for most of the dSphs. 
For example, ISchweitzerl jl996l) conclude that Ursa Minor 
appears consistent with it moving in the direction of the 
Magellanic C l ouds, although a more recent measurement by 
iPiatek et alJ J20O5) concludes that its membership of the 
Magellanic Stream is unlikely. Sculptor and Fornax have 
proper motion measureme nts consistent with their mem- 
bersh ip of the FLS stream iSchweitzerlll99(i iDinescu et alJ 
l2004f) , although an earlier measu rement for Fornax h ad ruled 
out membership o f this s tream JPiatek et alJl2002j) . 

lKroupa~ al. (2005) have suggested that the distri- 
bution of Galactic satellites is inconsistent with the dis- 
tribution of s atellit es expected from CDM simulations. 
IZentner"et al. (2005]), however, point out that a spherical 
distribution of satellites is not the correct null hypothe- 
sis for dealing with CDM simulations, and the results of 
iKroupa et alJ J2005:*) are marg inally consistent with a pro- 
late distribution of satellites. iLibeskind et alJ (2005) find 
that the subset of subhaloes which have the most massive 
progenitors at early times, and therefore arguably the ones 
most likely to be luminous, are distributed in a similar way 
to the Galactic satellites, even though the overall distribu- 
tion of dark m atter and dark sate llites is very different. In a 
similar way to lKnebe et al] (120041) . this result is interpreted 
in terms of the accretion of satellites along filaments. 



The satellite distribution of the Galaxy therefore 
plays an important role in cosmology and galactic 
dynamics. Whether this distribution is representative 
of satellite systems in general still needs to be an- 
swered. As a result, the analysis of the spatial distri- 
bution of similar systems is important to several fields. 
We are currently involved in an extensive stu dy of 
M31 and its environment dlbata et al] l200lt l2004 EoOE : 
Ferguson etafl 120021: IMcConnachie et al] l2003t l2004bic . 
20051: IMcConnachie fc Irwinl l2005bt IChapman et al] I200E : 
Irwin et al.ll2005l) . This galaxy hosts the next closest satel- 
lite subsystem to our own and in this contribution we con- 
duct an analysis of its spatial distribution and compare it 
to the Galactic satellites, to explore the uniqueness of our 
satellite system. Some of the res ults presented here have pre - 
vio usly been discussed briefly in | McConnachie et all (|2004a) 
and IMcConnachie fc Irwinl 12111 15a :. In Section 2 we discuss 
the initial set of galaxies which we study, incompleteness 
concerns, and the distance measurements that we use for 
each. We also define a M31 spherical coordinate system, in 
analogy to Galactic coordinates. In Section 3, we use the dis- 
tance and radial velocity data for the Local group galaxies 
to revisit the question of membership of the M31 subgroup. 
In Section 4 we analyse the spatial distribution of the satel- 
lite system and compare it with our own Milky Way. In 
Section 5, we search for "ghostly" streams around M31, in a 
similar way to lLvnden-Bell fc Lvnden-Belll lll995l) . Section 6 
discusses some of our findings, and Section 7 summarises. 



2 THE DATA 

Most Galactic satellites are located within some 200 kpc of 
the Sun. Distances to these objects are reliably determined, 
and are typically accurate to 5%. It is therefore relatively 
easy to conduct a reliable analysis of the distribution with 
respect to the Milky Way, as the data is of high quality and 
only a small correction needs to be applied to account for 
the offset of the Sun from the Galactic centre. 

M31 is located at a distance of ~ 780 kpc (eg. 
Freedman fc Madord ll99Ct iJoshi et al] 120031 : iBrown et all 
20041: IMcConnachie et aljl2005h . Some 14 to 18 other galax- 
ies are believed to be associated with it at roughly equivalent 
heliocentric distances (of order 600 — 1000 kpc). Accurate 
relative distances to each of these objects is crucial for a 
study of their distribution around M31, as the correction 
that must be applied to account for our offset from M31 is 
two orders of magnitude larger than for the Galactic centre. 
Additionally, it is important to have an understanding of any 
possible incompleteness in the dataset (ie. currently undis- 
covered group members) . In this section, we first present the 
data that we will be using, and define a coordinate system 
in which to analyse them. We then discuss the accuracy of 
the data, and comment on incompleteness effects. 

2.1 The Possible Members 

A full list of all galaxies which are initially considered as 
possible members of the M31 subgroup are listed in Table 1, 
along with their Galactic coordinates, adopted distance and 
uncertainty, and radial velocity. For completeness, we start 
by considering all Local Group candidates in the same part 
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Galaxy 




1 




b 


r (kpc) 


vq (bus 1 ) 


M31 


121. 


.2 


-21 


.5 


785 


± 


25 


-301 


M33 


133. 


.6 


-31 


.3 


809 


± 


24 


-180 


NGC 205 


120. 


.7 


-21 


.7 


824 


± 


27 


-244 


NGC 147 


119. 


.8 


-14 


.3 


675 


± 


27 


-193 


NGC 185 


120 


.8 


-14 


.5 


616 


± 


26 


-202 


And I 


121. 


.7 


-24 


.9 


745 


± 


24 


-380 


And II 


128. 


.9 


-29 


.2 


652 


± 


18 


-188 


And III 


119. 


.3 


-26 


.2 


749 


± 


24 


-355 


AndV 


126. 


.2 


-15 


.1 


774 


± 


28 


-403 


And VI 


106 


.0 


-36 


.3 


783 


± 


25 


-354 


And VII 


109. 


.5 


-9 


.9 


763 


± 


35 


-307 


And IX 


123 


.2 


-19 


.7 


765 


± 


24 


-216 


LGS3 


126 


.8 


-40 


.9 


769 


± 


23 


-286 


Pegasus 


94. 


.8 


-43 


.5 


919 


± 


30 


-182 


WLM 


75. 


.9 


-73 


.6 


932 


± 


33 


-116 


DDO210 


34 


.0 


-31 


.3 


1071 


± 


39 


-137 


M32 


121. 


.2 


-22 


.0 


785 


± 


25 


-205 


IC10 


119. 


.0 


-3 


.3 


825 


± 


50 


-344 


IC1613 


129 


.8 


-60 


.6 


700 


± 


35 


-234 


EGB0427+63 


144. 


.7 


-10 


.5 


1300 


± 


700 


-99 



Table 1. All candidate members of the M31 sub-group initially 
considered in this study, along with positional information and 
heliocentric radial velocities. M32 is assumed to be at the same 
distance as M31. The three objects listed below M32 have had 
their distances and uncertainties taken from Mateo (1998); all 
those objects listed above M32 have had their distances and un- 
certainties taken from Papers I & II. The heliocentric radial veloc- 
ities for all objects other than Andromeda IX are those listed by 
the NASA/IPAC Extragalactic Database (NED). Andromeda IX 
was only been recently discovered by Zucker et al. (2004) and 
its radial velocity has newly been measured by Chapman et al. 
(2005). 



Galaxy 


1 


b 


r (kpc) 


LMC 


280.5 


-32.9 


19 


SMC 


302.8 


-44.3 


58 


Sagittarius 


6.0 


-15.1 


24+2 


Ursa Minor 


105.0 


+44.8 


66+3 


Sculptor 


287.5 


-83.2 


79+4 


Draco 


86.4 


+34.7 


82+6 


Sextans 


243.5 


+42.3 


86+4 


Carina 


260.1 


-22.2 


101 ± 5 


Fornax 


237.1 


-65.7 


138 ± 8 


Leol 


226.0 


+49.1 


250 ± 30 


LeoII 


220.2 


+67.2 


205 ± 12 


Phoenix 


272.2 


-68.9 


445 ± 30 



Table 2. The satellites of the Milky Way and their Galactic 
coordinates. All distance information is taken from Mateo (1998). 

2.2 M31-centric coordinates 

It is convenient to use a M31-centric coordinate system 
for this investigation and we define a system that is di- 
rectly analogous to Galactic coordinates, but centered on 
M31. In Galactic coordinates, the position of M31 is given 
by (ijvf3i, &M3i, fJVf3i)- Another galaxy, 5*, has a position in 
Galactic coordinates of (l,b,r). The coordinates of S in a 
Cartesian coordinate system, centered on the Galaxy and 
aligned with Za/3i, are 

x = r cos (I — Im3i) cos b 
y = r sin (I- I M31 ) cost 

2 = rsinfo. (1) 

The coordinates of S in the M31-centric spherical coordinate 
system are (l',b',r'), and are given by 



of the sky as M31 that are not located closer to the Milky 
Way than t o M31. The majority of the d istance estimates are 
tak en from iMcConnachie et alJ (|2004cT) (hereafter Paper I) 
and lMcCormachie^t" alTil2005ir (hereafter Paper II) and will 
be discussed shortly. For those galaxies not analysed as part 
of our M31 survey (objects listed below M32 in Table 1), 
we take their distances from lMateol dl998l) . M32 is assumed 
to lie at the same distance as M31. The radial velocities 
for all objects other than Andromeda IX are those listed 
by the NASA/IPAC Extragalactic Database (NED). The 
ra dial velocity of Androm eda IX has been newly measured 
bv lChapman et all i2005l) as part of the M31 radial velocity 
survey using the Deep Imaging Multi-Object Spectrograph 
on Keck II. 

Table 2 lists the satellites of the Milky Way which will 
later be compared to the M31 sate llites. The distance infor- 
mation is taken from lMateol Jl998t) . Leo A is not considered 
a member of the Galactic satellite system due to its rela- 
tively large separation from t he Milky Way. The putative 
dwarf galaxy in Canis Major iMartin et al.ll2004) has not 
been included as its true nature is still the subject of some 
debate. The ultra-fai nt Galactic companion in Ursa Major 
JWillman et al]l2005l) is not included as it has only recently 
been discovered and lacks a precise distance measurement. 



/' = ill,.:! ( ^7 I + 180" 



r — 
where 



l 12 , 12 l2\ 

(x +y +z 



(2) 



x + r A/3 i sini 

y 

z' + r A/3 i cosi 



R„ (90 - i) Rj (90 - 9) 



Ry(-bMBi) yy J- ( 3 ) 

We take 6 = 39°8 to be the position angle of the semi-major 
axis of M31 measured east from north in Galactic coordi- 
nates, a nd i = 12°5 as the in clination of M31 to the line- 
of-sight <lde Vaucouleurslll958l) . R is the appropriate Carte- 
sian rotation matrix. The resulting coordinate system is di- 
rectly analogous to the Galactic coordinate system but is 
centered on M31: I 1 is a longitude measured around the disk 
of M31, where I' = is defined to be the longitude of the 
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Milky Way, and b' is a latitude, such that b' = corre- 
sponds to the plane of the disk of M31. r' is the distance 
of S from the centre of M31. The Galaxy therefore lies at 
(/' = 0,6' = -12°5,r' = 785 kpc). Henceforth, the prime no- 
tation is dropped and the context will distinguish Galactic 
and M31-centric longitude, latitude and distance. The un- 
certainty in the distance to each candidate satellite trans- 
lates into uncertainties shared by all three M31-centric coor- 
dinates, and whose relative sizes depend upon the location 
of the individual galaxy relative to M31. 

2.3 Distance measurements 

A study based on distances to multiple objects benefits from 
having small random errors associated with each measure- 
ment. It is then straightforward to account for the size of the 
error in the analysis procedure and assess the confidence of 
any results. However, for these methods to be reliable it is 
crucial that there are no systematic differences between the 
measurements which remain unaccounted for by the formal 
uncertainty. 

Differential systematic uncertainties arise due to the 
heterogeneous datasets that are employed, the standard can- 
dle utilised, the observing/data reduction strategy and the 
algorithms applied to the data. Distances in astronomy are 
historically unreliable, even for the relatively nearby galax- 
ies of the Local Group. It is important, therefore, to use 
a homogeneous set of distances which have differential sys- 
tematic effects minimised. Any systematic uncertainties that 
may be present will then, presumably, affect each distance 
measurement in a similar way. 

The I-band magnitude of the tip of the red giant branch 
(TRGB) is an observationally and theoretically well deter- 
mined standard candle for relatively old, metal poor stel- 
lar populations. Using the Isaac Newton Telescope Wide 
Field Camera (INT WFC), we have obtained Johnson V 
and Gunn i photometry for a large number of Local Group 
galaxies, the majority of which are members of the M31 sub- 
group. In Paper I, we developed an analysis procedure to 
minimise the observational difficulties in measuring the lo- 
cation of the TRGB in a photometric dataset. Paper II then 
went on to apply this procedure to 17 Local Group galax- 
ies for which we have INT WFC photometry and derived 
a homogeneous set of distances, as well as metallicities, for 
these objects. The identical data acquisition, reduction and 
analysis procedure for each galaxy goes a large way towards 
minimising differential systematic uncertainties, and the for- 
mal uncertainties on the measurements are ~ 5 %, similar to 
the Galactic satellites. The measurements are all consistent 
with previous estimates for each of the targeted systems, 
and there is no evidence for any systematic offsets or biases. 
The reader is referred to these papers for a more thorough 
discussion of these measurements. For the purposes of this 
study, there are a few systems that need to be considered 
for which we have not derived distances in this manner and 
it is necessary to supplement our data with that from other 
studies. 

The distance to M31 is the key distance in this study. 
The distance calculated for it in Paper II is 785L25 kpc. This 
measurement made use of stars in an annulus sampling the 
halo region of M31, minimising the influence of the young 
stellar populations and the substructure in this galaxy. The 



annulus was far enough from the centre of M31 so as not 
to be affected by crowding, significant extinction or serious 
contamination from disk and bulge components. Instead, it 
sampled a predominantly old, metal-poor population and 
the TRGB was obvious as an abrupt change in star counts 
in the luminosity function. Any line-of-sight effects through 
the halo would be expected to smear out the TRGB, which 
was not observed (if this effect was significant, then we would 
expect instead to have measured the TRGB for stars in the 
halo of M31 which are closer to us, and so the distance to 
M31 would be underestimated ). The good a greement with 
recent RRLy rae (794 ± 37kpc ; lBrown et al.ll2004) . Cepheid 
(791 ±40 kpc: IJoshi et al.l2003l) an d independent TRGB dis- 
tance measurements f783±43kpc: IPurrell et al.ll200l|) sug- 
gests that our distance estimate is robust. 

2.4 (In) completeness of the M31 subgroup 

The pioneering surv ey by Sydney van den Bergh 
llvan den Bergblll972 a b) found the first dSph companions 
to M31: Andromeda I, II and III. This photographic (II- 
Iaj) survey was designed to detect objects fainter than the 
limits of the Palomar Sky Survey, and was extended in 
1973 - 74 to encomp ass a total area of ~ 700 sq. degrees 
llvan den Bergbl Il974f) . No new dwarfs were discovered in 
this region. The limiting surface brightness to which this 
survey is complete is unknown, and the sky coverage be- 
comes patchy beyond a projected distance from M3 1 of 
~ 100 kpc. Subsequentlv.lArmandroff et alJ Il998lll999l) and 
iKarachentsev fc Karachentseval 11999 ) examined POSS II 
Illaj survey plates in the vicinity of M31 and discovered 
three new companions to M31 (Andromeda V, VI, VII). 
During the same period, IWhiting et alJ ll997l I1999F) con- 
ducted an all-sky search for new Local Group galaxies. While 
this search found two new Local Group galaxies, no new 
satellites of M31 were discovered. 

The recent surveys for M31 satellites have been system- 
atic in their sky coverage surrounding M31 and the discovery 
of bodies such as Andromeda V bv lArmandroff et alJ <ll998T) 
(with a central surface brig htness ~ 25.3 mags arcsec -2 ; 
iMcConnachie fc Irwinl2005bh show that they have been sen- 
sitive to low surface brightness features. It seems probable 
that all the dwarf galaxy companions to M31 with prop- 
erties similar to those discovered prior to 2004 have been 
discovered, due to the large (~ 1500 sq. degrees), contiguous 
nature of the surveys. Some satellites, however, may lie di- 
rectly behind the disk of M31. The area of sky covered by 
the disk of M31 behind which satellites could not be ob- 
served is ~ 5 — 10 sq. degrees; only one or two satellites are 
likely to hide from detection by this means. 

The discovery of Andromeda IX, with a central sur- 
face brightness in V of 26.8 mags sq. arcsec -2 IZucker et alJ 
l2004h illustrates that the faint end of the M31 satellite lu- 
minosity function is yet to be fully explored. It is currently 
impossible to say if Andromeda IX is unique, or whether it 
is the first of many extremely low surface brightness M31 
companions that await to be discovered. In the hierarchi- 
cal CDM framework, the luminous satellites are expected to 
make up o nly a small subset of the total number of dark mat- 
ter haloes (Kauffm ann et al"I IT993') , with baryonic infall and 
star formation being suppressed in the majority, perhaps by 
reionisation or supernovae feedback (eg. iBullock et aljEoOd : 
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Figure 1. The influence of the Galactic disk on the detection of M31 satellites, as seen in the M31-centric coordinate system. Contours 
indicate the volume of space obscured by the Galactic disk. Objects residing in the region of sky indicated by the central contour will lie 
in the Galactic plane < 15° in Galactic coordinates) if they are > 100 kpc from M31; objects lying inside the next contour will lie in 
the Galactic plane if they are > 150 kpc from M31, and so on out to > 300 kpc from M31 (outer contour). Shrinking of the contours by 
a factor of 2 — 3 provides a good representation of the region of sky covered by the innermost b = ±5° of the Galactic disk, where we 
expect the extinction to be most extreme. 



iDekel &Ts ilk 1986). Cosmological simulations of galaxy for- 
mation have difficulty modelling all these processes, and star 
formation is usually dealt with by semi-empirical techniques. 
It is important to constrain the relative importance of these 
effects by obtaining an observational census of the number of 
these very faint luminous satellites around the Galaxy and 
M31. Severa l searches using SPSS are currently unde rway to 
do just this dZucker et all2004lwillman et all200Sl) . While 
these surveys are not yet completed, the results so far indi- 
cate that there is not a significant extra population of these 
objects. 

To the north of M31 lies an area of sky obscured by 
the Galactic disk, which could hinder the detection of M31 
satellites. Figure Q shows an equal-area Aitoff projection of 
the satellites of M31 in the coordinates defined previously. 
Due to the orientation of M31 with respect to lines-of-sight 
through the Galactic disk, obscuration of M31 satellites is 
a function of position vector from M31. We consider the 
Galactic plane to defined by |6| < 15° in Galactic coordi- 
nates. Objects which reside in the region of sky indicated by 
the central contour will lie in the Galactic plane if they are 
> 100 kpc from M31; objects lying inside the next contour 
will lie in the Galactic plane if they are > 150 kpc from M31, 
and so on out to > 300 kpc from M31 (outer contour). 

Several galaxies reside in the most badly affected quad- 
rant and a few of these, including one of the most intrin- 
sically faint, lie in or near the Galactic plane in areas of 
significant extinction, namely AndV, And VII, NGC 185, 
NGC 147 and IC10. For the first four of these galaxies, the 
level of extinction is modest, at 0.4 — 0.7 mags in V. Only 
galaxies which were already intrinsically very faint (compa- 
rable to And IX) might be hidden by this level of extinction. 
On the other hand, IC10 is in a very heavily obscured re- 
gion of the sky and , according to the extinction maps of 
ISchleeel et alJ l)l998t) . has over 5 mags of extinction in V. 
The exact value is unreliable, given the patchy, irregular 
and large amount of extinction very close to the plane of the 



galaxy. If other M31 satellites lie in these heavily-obscured 
regions, it is not clear that they would have been easily ob- 
served; IC10 lends itself to discovery as it is currently under- 
going a starburst. If the contours in Figure are shrunk by 
a factor of 2 - 3, they provide a good representation of the 
area of sky covered by the innermost b — ± 5° of the Galac- 
tic plane. This is the region of the Galactic plane where we 
are most likely to get these extreme reddening values. Con- 
sidering the small area affected, it is unlikely that more than 
a few satellites (at most) will be affected. 



3 MEMBERSHIP OF THE M31 SUBSYSTEM 

A useful definition for the membership of the M31 satellite 
system are those objects whose kinetic energies are insuffi- 
cient to escape the gravitational potential of M31. In addi- 
tion, their separation must be small enough such that it is 
more likely for them to be bound to M31 rather than the 
Local Group as a whole. Choosing the correct membership 
is complicated by uncertainties in the mass and extent of the 
M31 halo, the separation of each candidate from M31 and 
the velocity of each candidate relative to M31. However, the 
uncertainty in the separation of the candidates from M31 is 
generally < ±40 kpc and is sufficient for this purpose. Like- 
wise, the mass of M31 , although not known precisely, is of or- 
der 1 - 2 x 10 12 Mrr. Evans fc WilkinsonllgOOri lEvans et alJ 
l200Ct llbata et alJl2004 iFardal et al.ll2005l) and is sufficient 
to ascertain membership. The greatest practical uncertainty 
in determining membership of the M31 subgroup comes from 
the velocity estimate. Only the heliocentric radial velocity 
of M31 and its companions can be measured directly and 
from this an estimate needs to be made of their relative 
3-di mensional velocity. 

lEinasto fc Lvnden-Belll l|l982l) demonstrated that the 
tangential velocity of M31 is small in comparison to the 
radial velocity under the assumption of no net angular mo- 
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Figure 2. The distance of each candidate M31 satellite from M31 
compared with an estimate of the potential required for the satel- 
lite to be virially bound to M31, as described in the text. Solid 
horizontal lines represent the uncertainty in the separation of each 
satellite from M31. The horizontal dot-dashed line represents the 
escape velocity from an isothermal M31 halo of circular velocity 
~ 250 kms -1 . The dotted lines correspond to the escape velocity 
for a point mass of 1, 2 and 3 X 10 12 Mq. The vertical dashed line 
is the distance of M31 from the barycentre of the Local Group 
assuming the mass of M31 is of order the mass of the Galaxy. 
Galaxies lying in the area of the diagram that is approximately 
delimited by these lines are most likely satellites of M31. 



mentum in the Local Group. Using this hypothesis, the to- 
tal Galactocentric velocity of M31, VM3i, is well approx- 
imated by the Galactocentric radial velocity component, 
v r> M3l = —123 kms -1 . We define the Galactocentric ra- 
dial velocity of S, located at a position r s from the Galactic 
centre, to be v r , a . The component of the velocity of S with 
respect to M31 in the direction r s is then 



Av 



l>M3ir s .rM31 



(4) 



where r s , rA/3i are unit vectors. Independent velocity in- 
formation for S in the remaining two orthogonal directions 
is unavailable. We therefore assume equipartition of kinetic 
energy and multiply (Av) 2 by a factor of three. The re- 
sulting quantity is an estimate of twice the specific kinetic 
energy of the galaxy, corresponding to an estimate of the 
gravitational potential required for the putative satellite to 
be virially bound to M31. In Figure [5] we plot this quantity 
against the separation from M31 for each candidate satellite. 
The horizontal solid lines indicate the uncertainty in the sep- 
aration of each candidate from M31. The horizontal dashed 
line corresponds to the escape velocity from an isothermal 
halo with a circular velocity of ~ 250 kms . The dotted 
lines correspond to the escape velocity for a point mass of 
1, 2 and 3 x 10 12 M . While these mass profiles are unrealis- 
tic, they represent simple models which bracket the extremes 
of the global potential, without being too prescriptive. The 
vertical line is the distance of M31 from the barycentre of 
the Local Group, assuming that the masses of the Galaxy 
and M31 are roughly equivalent. Objects lying in the area 
defined by these lines are probably satellites of M31. 



Figure H shows that M32, NGC 205, 147, 185, An- 
dromeda I, II, III, V, VI, VII and IX, LGS3, IC10 and M33 
are likely bound satellites of M31. The membership of Pe- 
gasus and IC1613 to the M31 subgroup is marginal as they 
are located ~ 500 kpc from M31. Nevertheless, we consider 
them to be members as their current position relative to M31 
and the Galaxy means that their orbits are almost certainly 
dominated by the potential of M31. WLM and Aquarius are 
located ~ 1 Mpc from M31, and should not be considered 
as satellites. The separation and velocity of EGB0427+63 
means that it, too, is unlikely to be a satellite. 



4 THE 3D DISTRIBUTION OF THE 
ANDROMEDA SUBGROUP 

Figure |3] shows the 3-dimensional distribution of the 16 
members of the M31 subgroup. The plane indicated is 
aligned with the disk of M31 and each grid cell is 100 x 
100 kpc. Solid lines are used for satellites that lie above the 
plane of the disk and dotted lines are used for satellites that 
lie below the plane of the disk. The direction of the Galaxy 
is marked. The top panel of Figure 0| shows the same distri- 
bution as an equal-area Aitoff projection of the M31-centric 
I and b coordinates. The position of the Galaxy is included 
as a reference. The lower panel shows error bars which repre- 
sent the combined uncertainty in the position of each satel- 
lite due to the uncertainty in its distance and in the distance 
to M31. The positions of M32, NGC 205 and AndromedalX 
in this projection are sensitive to these uncertainties due to 
their close proximity to M31, but the positions of the ma- 
jority of the satellites are relatively well constrained. 

Several general points are clear from inspection of the 
top panel of Figure 0| For example, half of the satellites lie 
within \b\ < 20° of the disk of M31. The satellite distribution 
as a whole does not appear isotropic; the upper left quad- 
rant and lower right quadrant house only one satellite each, 
while the remaining fourteen are confined to the other two 
quadrants. In addition, the satellites appear to congregate 
within 90° of I = 0. 



4.1 Variation with latitude 

In the lower panel of Figure^ the differential absolute (M31- 
centric) latitude distribution of the M31 satellites is plot- 
ted. In the top panel, we show the Galactic satellite dis- 
tribution as a comparison (in Galacti c coordinates). The 
putative dwarf galaxy in Canis Major iMartin et al.ll2004h 
has been highlighted in the upper panel as a dot-dashed 
line. The dashed curves correspond to an isotropic distri- 
bution of satellites. The correct null hypothesis for a CDM 
distribution of satellites may not be isotropy, however, as 
the satellites are ex pected to have a prolate distribution 
JZentner et alJl2005ft . This makes no difference to the fol- 
lowing discussion. 

With the inescapable caveat of small-number statistics, 
the M31 satellite distribution is qualitatively similar to the 
isotropic distribution and suggests that isotropy is a reason- 
able first-order assumption for the latitude distribution of 
satellites. Formally, a Kolmogorov-Smirnov (K-S) test be- 
tween the Galactic and M31 distributions shows that the 
two populations are not significantly statistically different. 
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Figure 3. A three-dimensional representation of the spatial distribution of the M31 subgroup. The plane shown is aligned with the disk 
of M31. Each grid cell is 100 X 100 kpc. I and b are the M31-centric spherical coordinates defined earlier. Solid lines are used for satellites 
lying above the disk, while dotted lines are used for satellites lying below the disk. The direction of the Galaxy is marked. 
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They are consistent with being drawn from the same under- 
lying distribution at the 38 % level. The inclusion of Canis 
Major makes the distributions consistent at > 60 % proba- 
bility. However, the Galactic system has an apparent under- 
abundance of satellites at low latitudes in comparison to 
the isotropic case and to M31; only one (definite) satellite is 
observed around the Galaxy at \b\ < 20° (Sagittarius) com- 
pared to seven around M31. If the Galactic satellites follow 
an isotropic distribution, then we would expect ~ 4 satellites 
to be located at \b\ < 20°, instead of the 1(2) so far discov- 
ered. This implies that obscuration by the Galactic disk at 
low Galactic latitudes may have led to an incompleteness 
in the Galactic subgroup of 15 — 25 %. While approximate , 
these values agree with the study bv lWillman et all i2004l) . 
who conclude that there may be a ~ 33 % incompleteness 
in the Galactic satellites caused by obscuration at low lati- 
tudes. 



O 1 1 1 1 1 1 1 1 1 

20 40 60 80 

M 

Figure 5. The differential absolute latitude distribution of the 
Galactic (upper panel) and M31 (lower panel) subgroups. The 
dashed curves shows an isotropic distribution for the same total 
number of satellites. Also shown as a dot-dashed line in the top 
panel is the position of the putative Galactic satellite in Canis 
Major. There is an apparent deficit of Galactic satellites at low 
latitudes, which is most readily explained as a selection effect 
due to obscuration by the Galactic disk. The M31 and Galactic 
latitude distributions are not statistically different however; a K- 
S test shows they are consistent with being drawn from the same 
underlying distribution at the 38 % level. 



4.2 Radial distribution 

The top panels of Figure [S] show the differential radial distri- 
bution of Galactic satellites (left) and M31 satellites (right). 
The lower panels show the respective cumulative distribu- 
tions. The abscissa have not been normalised by the virial 
radii of the host haloes because these are uncertain and are 
believed to be similar for the two hosts (25 8 kpc for the 
Galaxy, 280kpc for M31: iKlypin et alJl2002l) . The Galac- 
tic subgroup is noticeably more centrally concentrated than 
the M31 subgroup; half of the Galactic satellites are within 
~ 100 kpc while the corresponding separation for the M31 
system is ~ 200 kpc. However, the distributions are not 
statistically different; a K-S test on the distributions, nor- 
malised to the same total number of satellites, shows there is 
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Figure 4. Top panel: An Aitoff equal-area projection of the M31 satellite distribution, in the M31-centric spherical coordinates defined 
earlier. The position of the Galaxy is shown for reference and is defined to lie at I = 0. Bottom panel: Same as top, but now with error 
bars showing the combined uncertainty in the position of each satellite, due to the uncertainty in its distance and in the distance to M31. 
The positions of NGC205, M32 and Andromeda IX are relatively uncertain due to the close proximity of these objects to M31. All of 
the other satellites occupy relatively well-defined positions. 



a 24 % chance that they are drawn from the same underlying 
distribution. 

In a recent study, [Willman et alJ <l2004l) compared a 
simulated CDM radial distribution of satellites with those 
for the Galaxy and M31. These distributions were compared 
by assuming that dwarfs inhabit the sub-haloes with the 
highest circular velocity at zero redshift, and assuming that 
the Galactic satellite system is incomplete beyond 100 kpc. 
Based upon this comparison, they concluded that the Galac- 
tic system is potentially incomplete by a factor of up to 3 
(including the effects of incompleteness at low Galactic lati- 
tude), and that there may be some satellites beyond 100 kpc 
which have eluded detection, with absolute magnitudes and 
central surface brightnesses similar to the detected popula- 
tion. 

As the discov ery of the ultra-faint stellar system in Ursa 
Major illustrates jWillman et al.l200rf) . it is almost certainly 
the case that very faint satellites still await to be discovered 



around the Galaxy. It is also likely that the Galactic disk 
hinders the detection of low latitude Galactic satellites. It is 
not so clear, however, that there i s significant incomplete- 
ness in the outer (> 100 kpc) halo. Ilrwinl l)l994f) conducted 
an automated search of Schmidt sky survey plates for Galac- 
tic satellites, covering ~ 2/3 of the sky and extending down 
to a Galactic latitude of |6| = 20°. This survey was sensitive 
enough to detect systems such as Drac o and Ursa Mino r 
out to 300 kpc. Only Sextans was found l|lrwin et alJll99Ctl . 
The survey could have found systems up to one magnitude 
fainter than Sextans out to a distance of 200 kpc. Despite 
analysing more than 1000 sky survey plates, covering both 
hemispheres, no other comparable satellites were found. In 
addition, given that the total numbers of satellites around 
M31 and the Galaxy are roughly equivalent, and that by 
normalising the populations by total number the distribu- 
tions are broadly similar, it seems unlikely that the Galactic 



M31 7 s Satellite Galaxies 9 



Galactic subgroup 
12 members 



E ° 



M31 subgroup 
16 members 




100 200 300 400 
r (kpc) 




200 300 
r (kpc) 



Figure 6. The differential radial distribution of Galactic satellites (top left) and M31 satellites (top right). The lower panels show the 
respective cumulative distributions. The Galactic subgroup is more centrally concentrated than M31; roughly half of its members are 
located within 100 kpc compared to ~ 200 kpc for M31. A K-S test, where the populations are normalised to the same total number of 
satellites, shows the distributions are not significantly statistically different; formally, there is a 24% chance that they are drawn from 
the same underlying distribution. 



po pulation is seriously in complete in the manner suggested 
bv lWillman et alJ J2004h . 



4.3 Distribution with respect to the Galaxy 

Figure |7| is the Aitoff projection shown in Figure 2] with a 
dashed line representing the plane separating the near side 
of M31 from the far side, with respect to the Galaxy. Ob- 
jects lying inside this line lie on the near side of M31, while 
objects lying outside this line lie on the far side of M31. Ne- 
glecting distance uncertainties, NGC205 is the only object 
that clearly lies on the far side of M31; the other satellites 
are all close to, or in front of, the plane. It is important to 
emphasise that this plane is not a fit, and is defined only 
by the position of the Galaxy with respect to M31. In fact, 
fourteen satellites lie in front of the plane and two behind 
it. A suitable null hypothesis is that a satellite is equally 
likely to lie in front or behind this plane. This is easily tested 
with binomial statistics and is inconsistent with the observed 
distribution at 99.8 %. This does not take into account the 
distance uncertainty associated with each satellite. There- 
fore, we have also adopted a Monte-Carlo sampling tech- 
nique where we continually generate new M31 satellite dis- 
tributions by choosing the distance of M31 and its satellites 
from Gaussian distributions centered on the distance given 
for them in Table 1, with a standard deviation equal to 
the corresponding distance uncertainty. Only in 0.5 % of the 
generated cases do we find a distribution with equal num- 
bers of satellites in both hemispheres, indicating that the 
anisotropy is not an artifact of the distance uncertainties. 

Figure|5]shows the positions of the satellites in a Carte- 
sian coordinate system, such that the abscissa is aligned with 
the line-of-sight to M31. (xm31 = 0, j/M3i = 0, ZM31 = 0) is 



the position of M31 and (xm3i = 785, J/M31 = 0,za/3i = 0) 
is the position of the Galaxy. The distribution in the top- 
panel is affected most by the distance uncertainties, and the 
dot-dashed line represents the corresponding histogram for 
the Monte-Carlo distributions generated previously. 

For an isotropic distribution of satellites, the geometric 
centre of the satellite distribution is expected to be coin- 
cident with the position of the host. It is straightforward 
to calculate the mean (or median) and random error in 
the mean (a = a^/yN, where ad is the standard devi- 
ation) for the three distributions shown in Figure 8. The 
random error in the mean is the relevant statistic to use in 
this individual uncertainties in the position of each 

data-point statistically average out when looking at the dis- 
tribution of all data points (assuming these uncertainties are 
of comparable magnitude and are random). As the distribu- 
tions deviate from Gaussians, we have calculated ad using 
the more robust median- of-the-absolute-devia tion-from-the- 
median (MAD) statistic llHoaglin et al.llT98ot) . This is given 
by 

MAD = median \xi — median (x)\ (5) 

and is related to ad by ad = 1.48 x MAD for a Gaussian. 
This statistic gives an estimate of the scatter in a distri- 
bution even if it deviates significantly from a Gaussian and 
contains several outliers. We have also calculated ad using 
the mean-of-the-absolute-deviation-from-the-mean statistic, 
finding good agreement with the MAD statistic. 

The mean, median and a for all the M31 satellites are 
listed in Table [3] Also listed are the equivalent numbers 
when the outlying satellites Pegasus and IC1613 are ex- 
cluded. Table 2] lists the equivalent results for the Galac- 
tic satellite system, where the coordinate system has been 
shifted to be centered on the Galaxy. For this system, the 
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Figure 7. Same as Figure [I] but with a dashed line indicating the plane which separates the near side from the far side of M31. The 
satellites of M31 show a correlation with this plane, such that they preferentially lie on the near side of M31 with respect to the Galaxy. 
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Table 3. The mean, median, and random error in the mean (<x) for the three distributions shown in Figure|S] These have 
been calculated for the case where all the satellites are included and for the case where the outlying Pegasus and IC1613 
dwarfs are excluded. 



iLvnden-Bell fc Lvnden-Belll <ll995l) . This requires finding 
satellites that align on great circles, which may reflect a 
stream-like origin and assumes that there has not been sig- 
nificant precession of the orbits which would wash out any 
alignments, as would be expected if the halo of M31 is flat- 
tened. Here, we assume this is not the case. 

The probability of finding n satellites from a sample 
of N in one particular plane is readily estimated from bi- 
nomial statistics and is small if n is a large fraction of N. 
However, the probability of finding n satellites in any plane 
is much larger and reflects the fact that all possible planes 
will have been tried in the search for the plane around which 
the scatter is smalle st. This point is relevant to the claim by 
iKroupa et alJ J2005f) that the planar distribution of Galactic 
satellites is inconsistent with being drawn from a spherical 
distribution. This is only true at a high significance level if 
the plane in which the satellites align is considered to be 
predefined. 

Figure [5] is a Lambert equal-area zenithal projection of 
the polar paths of the satellites of M31. Each curve cor- 
responds to the track containing all points in the north- 
ern hemisphere perpendicular to the vector connecting the 
satellite to M31. Where each curve crosses another defines 
the pole of the plane containing the two satellites and M31, 
and is equivalent to the direction obtained on taking the 
cross-product of the position vector of the satellites from 
M31. Points at which multiple paths cross define the poles 
of planes which contain candidate streams of satellites. 

Several possible streams are highlighted by Figure[5]but 



random error on the mean is of a similar magnitude, or 
greater, than the mean/median of each coordinate. This sug- 
gests that there is no statistically significant offset between 
the centre of the Galactic satellite system and the Galaxy. 
Likewise, for the M31 subgroup, there is no statistically sig- 
nificant offset between the centre of its satellite system and 
M31 for yuz\ and zmh (taking into account the effects of 
the two outliers). However, the mean/median of xjv/31 is in- 
consistent with zero at the 3 a level. This suggests that the 
satellite distribution of M31 is significantly offset from M31 
in the direction of the Galaxy/barycentre of the Local Group 
at the 3<t level. 

It is pertinent to calculate whether an equivalent asym- 
metry exists in the velocities of the satellites. Using the no- 
tation defined in Section|H] the mean velocity of the satellites 
relative to M31 in the direction of the Galaxy is estimated 
as (Av) / (cos#M3i), where cos#m3i = r s .rM3i- This is the 
projection of the relative velocity of the satellite with respect 
to M31 onto the line-of-sight to M31, and is +19T21 kms" 1 . 
A positive sign indicates that the satellites are approaching 
more slowly than M31. The uncertainty is the random er- 
ror in the mean, and the velocity offset is not statistically 
significant. 

5 GHOSTLY STREAMS? 

Using the 3-dimensional information available on the M31 
satellite distribution it is possible to search for "ghostly 
streams" of satellites around M31 in a similar way to 
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Table 4. The same statistics as in TablelSlbut for the Galactic satel- 
lite system. The coordinates are the same as defined previously, but 
shifted to be centered on the Galaxy instead of M31. No statistically 
significant offsets are observed between the satellite distribution and 
the Galaxy. 




NGC205 



Figure 9. A Lambert zenithal equal-area projection. I increases anticlockwise as shown, b decreases radially from b = 90° at the centre 
to b = at the perimeter; the embedded circles represent 10° increments in b. Each curve represents the northern 'polar path' of the 
labelled M31 satellite ie. the track representing all points tangential to its direction from M31. The point at which two tracks cross defines 
the pole of the plane containing the two objects and M31. Points at which multiple tracks cross define the poles of possible streams 
containing the relevant satellites and M31. 



attributing significance to them without kinematic informa- 
tion is premature. For example, the projected position of 
M32 is sufficiently unreliable to neglect it in this analy- 
sis, and NGC 147 an d 185 are probably a binary system 
Jvan den BergbllTool . so should really be treated as a sin- 
gle point. 

However, there are numerous possible candidate 
streams indicated in Figure |UJ of which some of the most 
well defined are: 

• M33, Pegasus, IC10 (and possibly M32) 
(I ~ 180°, fo~ 12°); 

• Andromeda III, VII and LGS3 (Z ~ 140°, b ~ 23°); 

• NGC 185, 147 and 205 (and possibly Andromeda II 
(Z~74°,b~24°); 

• Andromeda I, IC10, 1613, NGC147 and 185 (and pos- 
sibly M32) (Z~80°,&~42°); 



• Andromeda II, IC10 and LGS3 (and possibly M32) 
(I ~ 50°, 6 ~ 30°); 

• Andromeda III, V and IX (I ~ 10°, b ~ 56°); 

• Andromeda V, VI and Pegasus (Z ~ 0, b~80°). An- 
dromeda IX and NGC 205 share a pole nearby at 
(Z~0,Z>~72°); 

• Andromeda I, II, III, and VI (Z ~ 352°, b ~ 38°) 

Without doubt many, if not all, of these associations 
are chance alignments; it is difficult to think of a plausi- 
ble physical association between M33, Pegasus, IC10 and 
M32, for example. However, a few streams could be viable, 
such as those defined by the dwarf ellipticals or some of the 
dwarf spheroidals. It is particularly interesting to note the 
planes defined by the latter group; Andromeda I, II, III, V, 
VI and IX all share poles in the same general vicinity of 
(I ~ 0, b ~ 60°) and the great circle that they loosely form 
is readily identified in Figure 2] These constitute six of the 
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Figure 8. Distribution of satellites along the line-of-sight 
to M31 (top panel) where {xm 31 = 0, Vm31 = 0, ZM31 = 0) 
is the position of M31 (dot-dashed vertical line) and 
( x M3l = 785, J/M31 = 0, 2^M3i = 0) is the position of the Galaxy. 
The satellite distributions in the two directions orthogonal to this 
(middle and bottom panels) are also shown. Distance uncertain- 
ties mostly affect the distribution in Xm31- To illustrate this, the 
blue dot-dashed curve shows the distribution of satellites along 
the line-of-sight to M31, derived by continually drawing the dis- 
tance of M31 and its satellites from Gaussian distributions cen- 
tered on their derived distance with a standard deviation equal 
to their distance uncertainty. 



seven dSph companions to M31 and have broadly similar 
star formation histories and colours (Paper II). A physical 
association between some or all of these objects may there- 
fore be plausible, and could represent analogous systems to 
the Magellanic or FLS streams proposed around the Galaxy. 

Further analysis of the planes listed above re- 
quires an investigation of the kinematics of the candi- 
date streams. The study of th e Galactic satellites by 
iLvnden-Bell fc Lvnden-Belll |l995) investigated the specific 
energy of the members of the proposed streams, under the 
assumption that they were physically associated with each 
other and had the same specific angular momentum, invari- 
ant with time. The derived specific energies of the putative 
members are then expected to be equal if the satellites are 
physically associated. The specific energy of each satellite is 
given by 



E= \v 2 G + \h 2 r' 2 -4> , (6) 

where vq is the Galactocentric radial velocity of the satel- 
lite, h is the specific angular momentum and <j> is the gravita- 
tional potential of the Galaxy at the location of the satellite. 
The first term represents the energy due to the radial mo- 
tion of the satellite; the second term is the rotational kinetic 
energy of the satellite in its orbit, where h has been assumed 
to be a common constant for all stream members. 

An equivalent kinematic analysis to that of 
ILvnden-Bell fc Lvnden-Belll 1119951) is not directly possible 
for the M31 satellites. The heliocentric radial velocities 
that are measured for these objects do not provide full 
information on the M31-equivalent of vg and therefore 
the first term of Equation HJ cannot be calculated. This 
will remain the case until proper motions for the M31 
satellites become available. The additional uncertainties 
on the correct form and magnitude of (j>, the distance 
uncertainties and the unknown contribution of interlopers 
to the make-up of the candidate streams (which may be 
as high as 100 %) also makes a numerical approach to 
this problem unsatisfactory without prior information on 
the orbits of at least some of the satellites. Without this 
information, any conclusions on the physical association, or 
otherwise, of the members of the M31 subgroup are likely 
to be highly vulnerable and uncertain. 



6 INTERPRETATION 

Using the distances derived in Papers I and II, we have found 
that the M31 satellites are distributed anisotropically along 
the line-of-sight to this galaxy, and are skewed such that 
most satellites are on the near-side of M31. We have quanti- 
fied this distribution and measure a > 40 kpc offset between 
the centre of the M31 satellite system and M31. This offset 
is in the direction of the Galaxy/Local Group barycentre, 
and is significant at the 3<r level. 

We have also examined the spatial distribution of the 
M31 subgroup using earl ier distance e stimates to M31 and 
its satellites tabulated in iMateol <199Sh (Andromeda V, VI, 
VII and IX are recent discoveries and are not included in 
this paper, so we kept their distances the same). The larger 
distance errors in most of the other measurements and the 
inhomogeneous nature of the distance estimates which make 
up this compilation, unsurprisingly dilute, although do not 
remove, the asymmetry observed in the distribution. 

One of the primary motives of this paper was to make 
use of more accurate and more consistent distance estimates 
to study the satellite distribution of the M31 system. These 
estimates are the key to the subsequent analyses (see Sec- 
tion 2.3) and a comparison of them with respect to pre- 
vious distance estimates is given in Paper II. The signif- 
icant benefit obtained by using this dataset is that all of 
the measurements are obtained on the same system and in 
this respect represent the most homogeneous Local Group 
distance estimates available. The distance estimate to M31 
is clearly the most critical and bringing M31 closer by 25 
- 60 kpc would remove the anisotropy. While this is a pos- 
sibility, it requires the TRGB measurement to M31 to be 
systematically offset relative to the other TRGB measure- 
ments. The TRGB analysis, however, was designed to min- 
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imise the probability of this occurring. The possible errors 
which could affect M3f in particular were discussed in de- 
tail earlier. I n addition, this inte rpretation would require 
other TRGB JPurrell et al.1 1200 ill. RR Lyrae ferown et alJ 
12004) and Cepheid llJoshi et al]l2003l) distances to M31 to 
be wrong in similar ways by similar amounts. If M3I is lo- 
cated at ~ 740 kpc, it is not obvious why all the recent 
distance estimates to it, which have made use of a range of 
independent standard candles, should produce very similar 
wrong answers. As such, we believe that this is the most un- 
likely interpretation of the result. It should not be neglected, 
but it potentially has much wider reaching consequences for 
the use of these distance indicators. 

Selection effects were discussed in Section 12.41 It is 
hard to understand how this result could be explained by 
this means. We are unaware of any effect which would 
prevent the detection of M31 satellites located between 
800 — 1000 kpc from the Galaxy, assuming that they plau- 
sibly have typical properties of the satellite system. These 
satellites would have to be located directly behind the disk 
of M31 to evade detection, and this is a very specific and 
small volume of space. Only ~ 5 — 10 sq. degrees of sky are 
masked by M31, compared to the ~ 1500 sq. degrees over 
which satellites have been surveyed for and discovered. 

It is an interesting question to ask if such an offset could 
have a dynamical origin. The main difficulty with a dynami- 
cal solution is timescale: the orbital (dynamical) timescale of 
a satellite at a distance of ~ 100 — 200 kpc from a ~ 10 12 Mq 
host is 1.6 - 4.5 Gyrs, meaning that satellites will have com- 
pleted > 3 orbits around M31 over the course of a Hubble 
time. Given this, it is unlikely that any anisotropy present in 
the satellite distribution could survive over a Hubble time, 
and suggests that it was either introduced recently or main- 
tained somehow. Further, any dynamical effect is presum- 
ably going to affect both the luminous satellites and the 
dark matter satellites, and so we would assume that the dark 
matter satellites must also be asymmetrically distributed in 
this manner. 

A scenario in which a satellite spends a larger fraction 
of its orbital period on the near side of M31 than the far side 
is easily envisaged in a simple point mass case. Here, if the 
major axis of the orbit of the satellite is aligned with the po- 
sition vector of the Galaxy, and M31 is at the far focus of the 
elliptical orbit, then the satellite is more likely to be observed 
on the near side of M31, as given by Kepler's laws. In reality, 
the potential of M31 is not that of a point mass and there is 
the added complication that the orbits of the satellites need 
to be correlated. If the M31 subsystem consists of dynami- 
cal groupings of satellites (Figure |3J|, then the probability of 
observing a heavily-skewed distribution might be expected 
to increase, as there are fewer independent satellite orbits. 
However, orbital precession in a non-spherical halo potential 
would destroy any initial correlations present. This may ei- 
ther indicate that many of the dwarf satellites of M31 were 
accreted relatively recently (within the past dynamical time 
or so) or that the outer potential of M31 is approximately 
spherical. There is a growing body of evidence from simu- 
lations of galaxy formation that the observed dwarfs con- 
sist of objects which were accreted relativel y recently (eg. 
iBullock fc Johnstonlliool lAbadi et alJl2005l) . else it is un- 
likely they could have survived tidal disruption until the 
present. This may help explain the skewed distribution. 



Another possible dynamical explanation could be that 
the M31 satellite galaxies are tracing the large-scale struc- 
ture of the nearby Universe. According to CDM structure 
formation, satellite galaxies are accreted preferentially along 
filaments. It may be, therefore, that the satellites were ac- 
creted along a filament that aligns with the Galaxy/Local 
Group barycentre, or that the satellites are tracing a dark 
filament that aligns with this direction. Alternatively, the 
potential around M31 may be asymmetric due to contribu- 
tions from other Local Group components e.g. a large Lo- 
cal Group halo and/or a contribution from the halo of the 
Galaxy, if it is significantly extended. This latter scenario 
will have been particularly relevant at early times when the 
separation of the M31 - Galaxy binary system was much less 
(the timing argument: iKahn fc Woltierlll959l : iLvnden-Belll 
1981) and if the Galaxy is significantly more massive than 
M31. However, effects such as precession could once again 
act to destroy any gross asymm etries present. 

Ivan den Bosch et al.l <|2005T) have recently shown from 
an analysis of the Two-Degree Field Galaxy Redshift Survey 
(2dFGRS) that, in general, the brightest galaxy in a group is 
offset in velocity from its satellites, and has a specific kinetic 
energy that is ~ 25 % of its satellites. For a relaxed CDM 
halo, this corresponds to an offset of the brightest galaxy 
from its satellite system of ~ 3 % of its virial radius, or 
of order 10 kpc for M31. The associated velocity offset is a 
few tens of kms . They suggest that the brightest galaxy 
either oscillates around the minimum of the relaxed halo 
potential (which is unlikely if the system has a cusp rather 
than a core), or that it resides at the minimum of the halo 
potential which is not yet relaxed. 

It is unclear whethe r the scenarios developed by 
Ivan den Bosch et alJ ll2005l) could apply to M31. These au- 
thors analysed groups of galaxies defined by four or more 
members in the 2dFG RS picked out usin g a group-finding al- 
gorithm developed bv lYang et al.l ll2005f) . The redshift range 
of these systems is 0.01 < z < 0.2. They found that for sys- 
tems where the velocity dispersion was less than 200 kms -1 
(corresponding to a total mass of < 1O 13 M0), the central 
galaxy was consistent with being located at rest at the cen- 
tre of the group. However, for this mass range, it was also 
consistent with being marginally displaced from the centre. 
It may be, therefore, that the anisotropic satellite distri- 
bution of M31 is tentative evidence that the halo of M31 
is not yet relaxed. To the best of our knowledge, there is 
currently no study which gives the virialisation timescale of 
halos formed in a CDM hierarchy as a function of mass. The 
rather unusual distribution of satellites of M31 that we have 
discovered surely merits further theoretical studies. 



7 SUMMARY 

We have re-derived membership of the M31 subgroup by 
examination of the positions and kinematics of the puta- 
tive members, and discussed their distance estimates and 
completeness. Using an M31-centric coordinate system, we 
have compared the distribution of satellites above and be- 
low the disks of M31 and the Galaxy. This suggests that the 
Galactic system is incomplete at low latitudes by ~ 20 %. 
Likewise, comparison of the radial distribution of satellites 
around these two hosts show that the M31 system is less cen- 
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trally concentrated than the Galactic system. We have also 
searched for "gh ostly streams" of satellites a r ound M31, in 
a similar way as iLvnden-Bell fc Lvnden-Bell l)l995h did for 
the Galactic system, and find several, including some which 
contain several of the dwarf spheroidal companions. Unfor- 
tunately, the current lack of proper motion data for these 
objects does not allow the possible physical association of 
these satellites to be better constrained. 

The analysis of the distribution of the M31 subgroup re- 
veals that the satellites are anisotropically distributed along 
our line of sight to M31, and the satellite system as a 
whole is skewed relative to its host in the direction of the 
Galaxy/Local Group barycentre. It is possible that this re- 
sult might be due to distance errors or selection effects, but 
such explanations would require unlikely fine-tuning and the 
former explanation in particular has potentially significant 
consequences for the reliability and use of the main distance 
indicators. Some possible dynamical explanations are sug- 
gested but present problems due to the timescales involved. 
It may be that many of the M31 satellites were accreted 
into M31's halo only relatively recently. Alternatively, it is 
possible that the asymmetric satellite distributi on that we 
observ e is related to similar findings bv lvan den Bosch et alJ 
(2005) at higher redshift, and that this may be evidence to 
suggest that the halo of M31 is not yet virialised. However, 
until an adequate explanation for this distribution is pre- 
sented, this finding warns against assuming that the M31 
subgroup represents a complete, relaxed and unbiased pop- 
ulation of galaxies. 
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